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Abstract For any positive integer n, let ¢(m) and S(n) be the Euler function and the Smaran- 


dache function respectively. In this paper, we use the properties and the curve figure of these 

two functions to study the solvability of the equation J` S(i) = o( 2a), and prove that 
i=1 

this equation has only two positive integer solutions n = 1,10. 
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81. Introduction and result 


For any positive integer n, the famous F.Smarandache function S(n) is defined as the 
smallest positive integer m such that n divides m!. That is, S(n) = min{m: m € N, n|m!}, 
where N denotes the set of all positive integers. From the definition of S(n), it is easy to see 
that if n = p{'p5?---p?* be the factorization of n into prime powers, then we have 


S(n) = qa res): 


It is clear that from this properties we can calculate the value of S(n), the first few values of 
S(n) are: S(1) = 1, S(2) = 2, S(3) = 3, S(4) = 4, S(5) = 5, S(6) = 3, S(7) = 7, S(8) = 4, 
S(9) = 6, S(10) = 5, ---. About the arithmetical properties of S(n), some authors had studied 
it, and obtained many interesting results. For example, Lu Yaming [2] studied the solvability 
of an equation involving the F.Smarandache function S(n), and proved that for any positive 
integer k > 2, the equation 

S(myz + m2 +--+ + my) = S(m1) + Sm) +--+ + Smg) 


has infinite group positive integer solutions (m1, M2, =, Mx). 
Jozsef Sandor [3] proved that for any positive integer k > 2, there exist infinite group 


positive integers (m1, M2, ++- , Mx) satisfying the inequality: 
Also, there exist infinite group positive integers (m1, M2, +- , Mk) such that 


S(mı + mo +--++ mz) < S(mz) + S(mz2) + +++ + S(me). 
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Rongji Chen [5] studied the solutions of an equation involving the F.Smarandache function 
S(n), and proved that for any fixed r € N with r > 3, the positive integer n is a solution of 


S(n)" + S(n)" T+- + Sn) an 


if and only if 
n = plp"! +p? +- +1), 


where p is an odd prime satisfying p”! + p"? +---+1|(p—1)!. 

Xiaoyan Li and Yanrong Xue [6] proved that for any positive integer k, the equation 
S(n)? + S(n) = kn has infinite positive integer solutions, and each solution n has the form 
n = pny, where p = kn, — 1 is a prime. 

For any positive integer n, the Euler function ¢(n) is defined as the number of all positive 
integers not exceeding n, which are relatively prime to n. It is clear that (n) is a multiplicative 
function. 

In this paper, we shall use the elementary method and compiler program to study the 
solvability of the equation: 


(1) +5) +---+5(n) = 0 (MED), (1) 


and give its all positive integer solutions. That is, we shall prove the following: 
Theorem. The equation 


S(1) + S(2)+---+S(n) =¢ (e4) 


2 


has and only has two positive integer solutions n = 1, 10. 


§2. Main lemmas 


In this section, we shall give two simple lemmas which are necessary in the proof of our 
Theorem. First we have the following: 
Lemma 1. For any positive integer n > 100, we have the inequality 


Z T’ n? 
ac TO 
= 5 S199 inn 
Proof. From the mean value formula of S(n) (See reference [7]) 
T? x? zr? 


we know that there exists one constant N > 0 such that 


n 2 2 2 2 2 2 
S sli) g T n | 1 T n 2 T n 
i=1 





S72 Inn 1199 12 Inn ` 11.99 Inn 


holds for all positive integer n > N. We can take N = 100 by calculation. This completes the 
proof of Lemma 1. 
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Lemma 2. For Euler function $(n), we have the estimate 








1 +1 1 
b (e ) = mn ) el. a aD 
ln (2 ln AE) 
Proof. Let n = pi'p3? pp” be the factorization of n into prime powers, then there 
always exist some primes p1, po, «+: ps such that pipo---ps >n. From [1] we have 
x 
S"inp=2+0 (=). 
log x 


psg 


by this estimate we know that 


lan < X inp; < 5 ln pi < ps < 2lnn. 
i=1 PiX<Ps 
Thus 1 1 
>-< 5 — <InInp, < lnln(2lnn). 
P Pi 


p|n Pi [Ps 





1 1 
Note that ¢(n) = nil (1 i, if mnt ) is even, then 
Pp 
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1 
If mei is odd, we can also get the same result. This completes the proof of Lemma 1. 


83. Proof of the theorem 


In this section, we shall complete the proof of our Theorem. First we study the tendency 
of the functional digraph 


K jaa 3 1 T g 
Goe ees n(m) 1199 ne 








By use of Mathematica compiler program we find that f(x) > 0, if z > 100754. 
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figure 1 


From the figure 1 we know that if n > 100754, then 


7 T? n? _n(n+1) s; 1 n(n + 1) 
i) < . near i 2 
D 11.99 mn 4 0m ran) i a( 2 ) 2) 





If x € (100754, +00), we use Mathematica compiler program to compute f (x), then we 
find that the derivative f'(x) is positive, so (2) is also true if x > 100754. 

Now we consider the solution of (1) for all n € [1, 100754]. By use of the computer 
programming language, we obtain that the equation (1) has no any other positive integer 
solutions except n = 1, n = 10. This completes the proof of Theorem. 


References 


1] Tom M. Apostol, Introduction to Analytic Number Theory, New York, Springer-Verlag, 
1976. 
2| Lu Yaming, On the solutions of an equation involving the Smarandache function, Sci- 
entia Magna, 2(2006), No. 1, 76-79. 

3] Jozsef Sandor, On certain inequalities involving the Smarandache function, Scientia 
Magna, 2(2006), No. 3, 78-80. 

4] Pan C. D. and Pan C. B., Elementary Number Theory, Beijing University Press, Beijing, 
1992. 








5] Rongji Chen, On the functional equation S(n)"+S(n)"~!+---+.$(n) = n, Smaraandache 
Notions Journal, 11(2000), No. 1-2-3, 128-130. 


Vol. 4 On the solvability of an equation involving the Smarandache function and Euler function 33 





6] Xiaoyan Li and Yanrong Xue, On an equation related to function S(n), Scientia Magna, 
4(2008), No. 1, 148-151. 

7] Y. X. Wang, On the Smarandache function, Research on Smarandache problems in 
number theory, II (2005), 103-106. 

8] Zhang Wenpeng, The elementary number theory, Shaanxi Normal University Press, 
Xi’an, 2007. 





The computing programme is given as follows if n € [1, 100754]. 
# include “stdio.h” 
# include “math.h” 
# define N 100754 
int S(int n) 
{int ret=1,num=n; 
unsigned long int nn=1; 
for(ret=1;ret<=n;ret++) { nn=nn«ret; 
if(an%num==0) break;} if (ret>n) ret=n; 
return ret; } 
int SumS(int n) 
{int ret=0,i; 
for (i=1l;i<=n;i++) ret+=S(i); 
return ret; } 
int coprime(int i,int n) 
{ int a=n,b=i; 
while(a!=b) { if(a==0) return b; 
if(b==0) return a; 
if(a>b) a=a%b; 


else 

b=b%a;} 

return a; } 

int Euler(int n) 

{int ret=1,i; 

for(i=2;i<n;i++) {if(coprime(i,n)==1) ret++;} return ret;} 
main() 

{ int kk; 


for(kk=1;kk<=N;kk++) if(SumS$(kk) ==Euler((kk*(kk+1) /2))) 
printf( “rusult is % d\n” ,kk); 
getch (); } 
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Abstract For orthogonal projective matrix R, i.e., R? = R and RT = R, we say that A is 
generalized Hermitian matrix, if RAR = A*. In this paper, we investigate the least residual 
problem ||AX — B|| = min with given X, B, and associated optimal approximation problem 
in the generalized Hermitian matrix set. The general expressions of the solutions are derived 
by matrix decomposition. 
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81. Introduction 


Some symbols and notations: Let C?"*” be the set of all m x n complex matrices with rank 

r, HC"*” be the set of all n x n Hermitian matrices. Denoted by At, A*, rank(A) the Moore- 
Penrose generalized inverse, conjugate transpose, rank of matrix A, respectively. Moreover, In 
represents identity matrix of order n, and J = (en, @n-1,°+: ,€1), e; E€ C” is the ith column of 
In. ||. || stands for the Frobenius norm. Matrix R € C”*” is said to be projective (orthogonal 
projective) matrix, if R? = R (R? = R and R* = R). 

Definition 1.1. If A € C”X”, we say that A is centro-symmetric matrix, if JAJ = A. 

The centro-symmetric matrix has important and practical applications in information the- 

ory, linear system theory and numerical analysis (see [1-2]). As the extension of the centro- 
symmetric matrix, we define the following conception. 

Definition 1.2. For given orthogonal projective matrix R € C?*", we say that A € C”*” 
is generalized Hermitian matrix, if RAR = A*. Denote the set of all generalized Hermitian 
matrices by GHC"*"”. 

In this paper, we discuss two problems as follows: 

Problem I.(Procrustes Problem): Given orthogonal projective matrix R € R”*”, and 

X,Bec"*™, find A € GHC”*” such that 
|| AX — B|| = min. 

Problem II.(Optimal Approximation Problem): Given M e C”*”, find A € Sp such 

that 
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|| M- All = min ||M — All, 
AESE 

where Spg is the solution set of Problem I. 

Obviously, when M = 0, Problem II is changed into finding the least Frobenius norm 
solution of Problem I. 

Many important results have been achieved about the above problems with different ma- 
trix sets, such as centro-symmetric matrixl3l, symmetric matrix!4—®!, R-symmetric matrix!®—7] 
and (R,S)-symmetric matrix!®] set. In this paper, we investigate the above problems in the 


generalized Hermitian matrix set by matrix decomposition. 


§2. Preliminary knowledge 


In this section, we discuss the properties and structures of (orthogonal) projective matrices 
R € C?*” and A € GHC™*". 
Denote s = rank(I — R), we know that r+ s = n since R? = R. Suppose that p1, po,..., Pr 
and q1,q2,---;@s are the normal orthogonal basis for range R(R) and null space N(R) of R, 
respectively. Let P = (p1, P2, ..., Pr) € C?*" and Q = (q1, q2,- --, qs) € C?”**, then 
P*P = l., Q*Q = Is, (1) 
RP = P, RQ =0. (2) 
Lemma 2.1.(see [9]) Let matrix A € C”*™ and its full-rank factorization A = FG, where 
Fecr*",Gecrx™ then A is projective matrix if and only if GF = Ip. 


Lemma 2.2. R €C?*” is projective matrix, then 


I 0 P 
R=(P Q) ae (3) 
0 0 Q 
-1 P 
where matrix (P Q) is invertible, and ( P Q) D ae 
Q 
If R is orthogonal projective matrix, we have 
I 0 P* 
R=(P Q) (4) 
0 0 Q* 


where {P Q) is unitary matrix. 


Proof. Assume that the full-rank factorization of R is R = PP, we obtain from Lemma 
2.1 and (1) that 
P=P*R, PP = L. (5) 


Similarly, if the full-rank factorization of I — R is I — R = QQ, we generate 
Q = QU- R), QQ = Ls, (6) 
since (I — R)? = I — R. Connecting with (1)(2)(5) and (6), we know that (3) holds. The 
equality (4) is obvious since R* = R. 
Lemma 2.3. Given matrices R as in (4) and A € GHC”*", then 
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G 0\ [P 
A= (P Q) , Y GeHC"™". (7) 
0 0} (Q 


Proof. According to Lemma 2.2 and Definition 2.1, it is clear that (7) holds. 
Lemma 2.3 indicates that arbitrary matrix M € C”*” can be written as 


Mı Mə P* 


ies a) M3 Mı Q* 


§3. The solutions of Problem I and II 


Given matrices X, B € C"*™, partition 
P* x P* Bı 
X= and B= l (8) 
Q* X2 Q* B2 
where X1, Bı € C”*™ and X2, By E C®*™, 
We need the following two lemmas derived from References [7] and [8], respectively. 
Lemma 3.1. Suppose that matrices X1, Bı in (8), then matrix equation A;X1 = Bı is 
consistent for A; € HCO'™*", if and only if BiXtX = Bı and X{B, = BX), the general 
solution is 
Ay = Â +, — XX} )K (U, — Xi XP), 
where Ay = (i. = 221) Bx} 4 (BX) (= 222), VK € HO. 
Lemma 3.2. Given matrices X4, Bı in (8), then 


min || GX; — Bı ||=|| Bid, - XF X) | 
Gecrxr 


if and only if G = B, Xj + Ko(I, — X, XF), Y Ko € C"*". 


According to Lemmas 3.1 and 3.2, we obtain 





Lemma 3.3. For the above given matrices X1, Bj, 


min || A,X, — Bı ||=|| Bid, — X7 X) | 
Ai€HC™xr 
if and only if 


XE BXY = XP XB XıXŤ, (9) 
and the expression of A; is the same as that in Lemma 3.1. 
Proof. | A, Xı = Bı |? =|| Bı = ByX{X + Bı X} Xı = A, Xı \|? 


=|| Bir — XT X1) |? + || Bixf X: - AX ||? 


Hence, the least residual can be attained only if By X T Xı = A,X, which is consistent for 
A; € HC™*" under condition (9) by Lemma 3.3. The proof is completed. 
Based on the previous analysis, Problem I can be solved in the following Theorem. 
Theorem 3.1. Given matrix R as in (4), X, B € C”*™ and the partition (8), then 
ga lA B|? = ||Bi(, — XT X)|? + || Ball”, (10) 
if and only if (9) holds, at this time 
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0) Ĝ + (I. —- XXPK, a) 0) [P 
0 0) \Q 


+ 
LX) YK € HO", 


A=(P (11) 


a + * 
where G = (I, — ©) B, x} + (Bı Xİ) (r - 

Proof. According to the unitary invariance of Frobenius norm, formulas (4) and (7), we 
obtain 


AX — Bl? 

2 
G o\ [P 

= (P Q) X-B 
Oo Oo} Lor 


E G 0 Xı Bı 
0 0 Xə Bo 
= ||GX1 — By|| ? + || Bell’. 


Therefore, the problem (10) is equivalent to the following least residual problem 
|| GX, — Bı ||. 


bo 











min 
From Lemma 3.3, we know that the Sac can be attained if and only if (9), and 
G= Ğ+ (I. =X) Kl — XıXŤ), 
where K € HC”*” is arbitrary. Submitting G into (7), then (11) holds. 

The following lemma stated from [6]. 

Lemma 3.4. Let L € C9%*™, A € 09%41, T € C™*™, and A? = A = A*, I? =T =I", 
then || L — ALT ||= Eo || L — ANT || if and only if A(L — N)T = 0. 

Let Sp be the solution set of Problem I. We can easily verify from its definition that Sz is a 
closed convex subsets in matrix space C”*” under Frobenius norm. The optimal approximation 
theorem!!®l reveals that Problem II has unique solution, which can be expressed in the next 
theorem. 


Theorem 3.2. Suppose that the given matrix in Problem II is 


M, M\ [P* 
M= (P Q) E Oa 
M3 Ma} \Q* 

then 

min || M—A|| (12) 

AESE 
if and only if 

Cth XXH AMi,- XXT) o) [P 
A=(P Q) ( 1X1) 1X7) . (13) 


0 0 Q“ 
where Ĝ is the same as that in Theorem 3.1. 
Proof. By using the unitary invariance of Frobenius norm and Theorem 3.1, we obtain 
M, M CLUS Rd =A). 0 
jm- ajp [1 ae) _ (@+G-- axe — xxt) 
M3; Ma 0 0 
= || (Mı - G) - (I, — XA XP) Kp — Xi XP) |? 
+ || Me ||? + || Ms |? + || Ma |, 
then the problem (12) equals to solve the minimum problem 
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j Mı — Ĝ) — (I; — XIX) ae een X, X”) ||. 
r&n, | hn — G) = ( 1X7) K( 1X7) | 


Moreover, since || M|? = || 





MiMi I? + || Homi \|?, hence the above minimum problem can be 


transformed equivalently as 





M +M* > 
min || (— 1 — G1) — (L, — ae a = A 
KEHO *Yr 2 
We further deduce from Lemma 3.4 that 
(O, —- XR - X Xf) = (I = XXt) xe), (14) 


submitting (14) into (11), we obtain (13). 
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